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Abstract 



The inverse problem which arises in the Camassa-Holm equation is revisited 
for the class of discrete densities. The method of solution relies on the use of 
orthogonal polynomials. The explicit formulas are obtained directly from the 
f^ ■ analysis on the real axis without any additional transformation to a "string" 

type boundary value problem known from prior works. 
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t^ ■ 1 Introduction 

O. 

The purpose of this paper and its sequel is to present an alternative solution to two 
inverse problems v^rhich appear in the context of Camassa-Holm type equations. We 
concentrate in this paper on the Camassa-Holm equation (CH)[3]: 

Ut - Uxxt + 3uUrc = 2UxUxx + UUxxx, (1-1) 

while in the sequel we study the inverse problem for another equation, discovered 
by V. Novikov (VN) [7]: 

Ut — Uxxt + 4li Ux = SuUxUxx + U Uxxx- (1-2) 

Both these equations admit another, often used, representation 

rrit + TTixU + 2mUx = 0, m = u — Uxx, (1-3) 



'Department of Mathematics, University of Isfahan, Isfahan, 81746-73441, Iran; 
k.mohajer@sci.ui.ac.ir 

^Department of Mathematics and Statistics, University of Saskatchewan, 106 Wiggins Road, 
Saskatoon, Saskatchewan, S7N 5E6, Canada; szmigiel@math.usask.ca 



for the CH equation, and 

mt + {iTixU + 'imux) u = 0, m = u — Uxx, (1-4) 

for the VN equation. The solution to the inverse problem presented in this paper 
is obtained directly on the real axis, instead of transforming the problem to the 
inhomogeneous string problem as was done in [T] and adapting the method used by 
T. Stieltjes in [8]. 

We give a self-contained presentation of the inverse problem for Equation ()2.6p . 
emphasizing the role of orthogonal polynomials (Proposition [TT]) and culminating 
in the elegant inverse formulas stated in Corollary [TH In the sequel, for the case 
of the VN equation (jl.2p . we will show that a different class of polynomials is 
germane to the inverse problem, namely a class of biorthogonal polynomials (Cauchy 
biorthogonal polynomials) studied in [2]. 

2 CH inspired inverse problem 

The well known Lax pair system for the CH equation can be written as follows (see 

{l - Dl - zm)ij = 0, 

(A + {u + l/z)Dx - Ux/2 - 1/z)^ = 0. 

We would like to point out that this Lax pair yields a slightly different normalization 
than in (II. 3j) . here, m = 2u — \uxx- If we set m{x,t) = X^"=i mj{t)5xj where mj{t) 
and Xj{t) are C°° functions, we obtain what is called the n-peakon Lax pair for the 
CH equation. We assume that for every fixed t > the function ip{x, t) is continuous 
on M and xi(0) < 0:2(0) < • • • < Xn{0). Now the first equation of the system (j2.ip 
implies that on the interval {xj,Xj^i), we have 

tl; = Aje"" + Bje-"" , (2.2) 

It is well known that requiring ip — )■ as \x\ —^ cxd is consistent with the system 
(j2.ip . Let's therefore assume that Aq = 1 and Bq = 0. Then from the continuity of 
■0 and the jump conditions imposed by the first equation of (|2.ip we have 

A,\ .. . . A 
'j 



s: =^^(^) ' '^^^^ (2-^) 



where 



and 



Clearly, 



K,{z) = T,{z)T,_,{z)...T,{z), (2.4) 



n{z)={l_iZt I'^'^T'^), l<k<n. (2.5) 



mfce^'^'^ 1 + f mfc / ' 
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Theorem 1. The spectrum of the boundary value problem 

{1 - dI - zm)i) = 0, V ^ 0, /or |x| ^ oo (2.6) 

is given by roots of An{z) = 0. 

In order to get more information about the spectrum we can rewrite equation 
as an integral equation 



ip{x,z) = z K{x,y)Tp{y,z)dM{x), 

where M is the distribution function of the measure m(x) and K{x,y) = 2^~ 

is the unique Greens function of 1 — D^ vanishing at ±00. One can easily write the 

explicit form of this integral equation as: 

n ^ 

tP{x,z) =z^-e-\^-''^\'>P{xj,z)mj, (2.7) 

i=i 

where, by assumption, all mj > 0. Moreover, evaluating at Xj, and introducing the 
matrix notation 

^ = {i;ixi,z),^Pix2,z),...,i^{xr„z)f, E= [^e-l^-^'il], (2.8) 

and 

P = diag(mi,m2,... ,m„), (2.9) 

we obtain equation (j2.7p in the form of the matrix eigenvalue problem: 

ip = zEPip. 

E is an example of a single-pair matrix (see [5]) which is oscillatory and so is EP. 
Oscillatory matrices form a subset of totally nonnegative matrices, which can be 
characterized as matrices whose all minors of arbitrary size are nonnegative. Oscila- 
tory matrices are special in that they have simple, strictly positive, spectrum [5]. 
We therefore have the following description of the spectrum of the boundary value 
problem (j2.6p . 

Lemma 2. An{z) = det(/ - zEP) = TT (l - — ) where < Ai < A2 < • • • < A„ 

Ai. 

l<fc<n 

We need one more fact. 

Theorem 3. The Weyl function W{z) = — XTIT ^'^ ^ (shifted) Stieltjes transform 
)f a measure dfi{x) = J21=i (^j^Xji '^j > ^> ^^0.^ ^•s 



o 



W{z) = j+ I -^dfiix) 
' z — X 



where 7 = / -dfi{x). 



Proof. This theorem follows, after a slight adjustment of notation, from Theorem 
4.2 proven in [1] and the observation that VF(0) = which is a consequence of 

Tk{o) = 1. n 

The continued fraction expansion ([B]) turns out to be at the heart of the inverse 
problem, in a full analogy with the inverse problem for an inhomogeneous string 
which was studied in the 1950's by M.G. Krein ([B], [1]). 

Theorem 4. 



A, 



« - _e2-n + 



-|mjie 2^" 



o^Xn pZXji — 1 _l_ 



-|m„_ie 2^'"-! + 



_|^,e-2., + 



Proof. Since 



we have 



(2.10) 



Bj ^"U::;,'' ^'-''^ 



-An\ m ( ^n-1 



Bn f m„e^^"^,,_i + (1 + f m„)B 



n-l 



A„ (1 - |m„)^„_i - |m„e 22:"B„_i' 
This equation is equivalent to 



(2.12) 



4^ = -e^-" + i . (2.13) 

^n 1 

— §m„e~2^" H 



2x , ^'^-i 



Hence, the proof is complete by induction. D 

Suppose that Pj/Qj (0 < j < 2n) is the jth convergent of the negative of 
the continued fraction we obtained above so that we are approximating the Weyl 
function W{z). Then we observe that 



Pi^ _ l + fm^ 



(2.14) 



and 

Q^ ~ 1- fm„(l -62(^-^-1-=^")) 
Therefore, it can be verified tliat 



P2 = -(e^""-e^""-i)Pi+Po, 



Q3 = -m„_ie-2^"-iQ2 + Qi, 



z 






P3= 2"^n-ie-^""-iP2 + A. 

These equations guide us to the following theorem 



Theorem 5. 






P2k = -(e2-"-^+i - e2-«-'=)P2fc-i + ^2, 





Q2k+1 ~- 


Z 

~ 2 




P2k+l = 


z 
" 2 ' 


Proof. By induction. 






Corollary 6. 







-'iXn-k 



Q2k + Q2k-l, 



mn-ke ^''"-''P2k+P2k-1- 



deg(P2fc+i) = deg{Q2k+i) = k + 1, 
deg(P2fc) = deg((52fc) = A;. 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



D 



PjQj^i - Pj^iQj = i-iy'\ l<j< 2n. (2.20) 

Proof. Use induction and the previous theorem. D 

Corollary 7. 



(2.21) 



In preparation for the inverse problem we introduce the following notation. 
Notation: Given a polynomial p{z) we denote by p[j] the coefficient at the jth 
power of z. Moreover, n — A; is abbreviated to k' . 

With this notation in place we have 

Theorem 8. 

Q2fc[0] = l, P2k[0] = e^^^' , 

k 

Q2k+i[o] = 0, P2k+i[o] = 1, ^2fc+i[i] = 1^2 "^r 

Proof. By induction. D 



The Weyl function of the n-peakon problem reads: 
Proposition 9. 

The inverse problem hinges on the following approximation formulas: 
Theorem 10. 



P2k{z) _ ^^ 1 
Q2kiz) Z^ 



^(•2) - TT^ = ^(^IfcTl)' ^ = 0,...,n-l, z^oo. 



Proof. By the previous theorem we have 

P2k+i _ P2k ^ 1 ^ n(_}__\ ^ 

Q2k+1 Q2k Q2k+lQ2k z^^+'^ 

Also, 

P2k P2k-l -1 /n,' -*■ ^ 

Q2fc Q2k-l Q2kQ2k-l ^ 

From these and the definition of the Weyl function we get 

Ti^/ \ P2k-1 P2n{z) P2k~l ^^/ 1 n 

Q2fc-1 Q2n(^) Q2fc-1 ^z2fc-l^' 



Similarly, we see that 



Q2k ^Z^ 



Wiz)-^ = 0{-^), z^^. 



(2.22) 



D 



Since Q2k+ii0) = we set Q2k+iiz) = zq2k+i{z), deg q2k+i = k. 
Proposition 11. (orthogonality) 

x^Q2kix)dfi{x) = 0, j = 0,...,k-l, k = 0,...,n-l (2.23a) 

f q2k+iix)dfiix) = 1, /c = 0,...,n-l, (2.23b) 

x-^ Q2k+i{x)dfi{x) = 0, j = l,...,k — l, k = 0,...,n — l. (2.23c) 

Proof. By the first approximation formula of Theorem [10] we have 

W{z)Q2k+i(z)-P2k+i{z) = 0{^), k = 0,...,n-l. (2.24) 



Therefore, 

z\W{z)Q2k+i{z) - P2k+i{z)) = 0{^), j = 0,...,k-2. 

z" J 

Since W{z) is analytic at z = oo so is the remainder and its order there is 0(-j). We 
can therefore use the residue calculus. To this end we choose a circle F whose radius 
is large enough to contain the support of ji defined in Theorem [31 An elementary 
contour integration implies then 

^ /J ^%^dM-)^. = 0, j = 0,...,fc-2. 
Consequently, applying Fubini's theorem and Cauchy's residue theorem we have 

x^Q2k+iix)dfi{x) = 0, j = 0,... ,k-2, k = 0,...,n-l. 
Similarly, from the second approximation formula of Theorem [TU] we have 

W{z)Q2kiz)-P2kiz) = Oi^), k = 0,...,n-l. (2.25) 

Thus, 

x-^Q2k{x)dfi{x) =0, j = 0, . . . ,k — 1, k = 0, . . . ,n — 1. 

Also we have 

Em^h±M_^h±M = o{^), k = 0,...,n-l. (2.26) 

z z z^+^ 

Since <52fc+i[0] = and P2fe+i[0] = 1, we have 

1 [ f Z~'^Q2k+l{z) 

' ' a^[x)az — 1 = 0. 



27ri ./r ./ z — x 

Hence, 

/ x~'^Q2k+i{x)dn{x) = 1, A: = 0, ... ,n- 1. 

Finally, using Q2k+i{x) = xQ2k+i{x) we obtain the claim. D 

Proposition 12. 

D . X f zQ2k{z) - xQ2k{x) fnn^ \ 

P2k{z) = / z r dfi{x), (2.27a) 

J x{z - X) 

r, / ^ [ zQ2k+l{z) - xQ2k+l{x) , . /n 0'7K\ 

P2k+i{z) = / ^ d^x{x). (2.27b) 

J x[z-x) 



Proof. Prom the second approximation formula of theorem [10] we have 

Wiz)Q2kiz) = P2kiz) + 0{^), z^oo. 
By Theorem El VK(z) = 7 + / -^-dfj,{x) which imphes 

D / \ n / \ i f Q2k{z) - Q2k{x) , ( . 
J z-x 

Recalhng that 7 = / -dfi{x) and performing elementary manipulations we obtain 
the claim. The proof for P2fc+i is analogous. D 

Let Cj = J x^diJ,{x), j G Z he the jth moment of the measure dfi{x) and let 
/ = [lij] = [cj+j] be the Hankel matrix of the moments of the measure dfi{x) 
appearing in Theorem [3l We note that the shift 7 = c_i. In addition, let us denote 
by A^ the minor of the submatrix of / of size kxk which starts at /oi- By convention, 

Ai = 1, J e z. 

Theorem 13. Suppose we are given an arbitrary discrete measure diJ,{x) = 
^ • aj6x,, Xj > 0, and define its shifted Stieltjes transform 

W{z)='y+ / dfi{x), 1=1 -d^{x). 

J z — X J X 

Let k be a natural number such that < k < n — 1. Then there exists a unique solu- 
tion to the approximation problem stated in Theorem [721 whose polynomials satisfy 
degree requirements dictated by Corollary^' 



Q2k{z) 



Ai 



1 z z ... z'^ 

Co Ci C2 ... Cfc 

Cl C2 C3 ... Cfc+1 

Cfc-1 Ck Cfc+l . . . C2k~l 



(2.28a) 



Q2k+1 



AO 



z z 

Cl C2 



z 

C3 



Ck Cfc+1 Cfc+2 



^fc+1 
Cfc+1 



C2k 



(2.28b) 



P2k{z) 



zQ2k{z) -xQ2k{x) 

x{z — x) 



dfi{x) 



p . X , zQ2k+i{z) - xQ2k+i{x) , . 
P2k+i{z) = / z ^ dn{x) 

X[Z — X) 



(2.28c) 
(2.28d) 



Proof. The proof is based on the work of Stieltjes [8]. We give only a complete argu- 
ment for Q2k, the proof for Q2k+iiz) is analogous. Both proofs rely on Proposition 

HI 

We write Q2k{z) = Yli=o1i^^- ^o^i using orthogonality conditions and the fact 
that Q2k[0] = 1 we have 

x^{l+^qix')dn{x)=0, j = 0,l,...,k-l, k = 0,...,n-l. 
-' i=i 

Recall that Cj = J x^dfi{x). Then we have 

k 

'^Ci+jqi = -Cj, j = 0,l,...,k-l. 



i=l 



Hence we obtain the system 
where 



Bq = -c, 



B 



(ci C2 
C2 C3 



Cfc \ 
Cfc+1 



{qi,q2---,qkV, c= (cq, ci, . . . ,Cfc_i)-' 



\Cfc Cfc+l . . . C2k-\j 

Thus, Cramer's rule implies that 



Q2k{z) 



Al 



1 z z^ 

Co Cl C2 
Cl C2 C3 



z 

Ck 
Cfc+l 



D 



Cfc-l Cfc Cfc+1 . . . C2fc_l 

The formulas for P2fc and P2k+i were already obtained in Proposition [T2j 

Using Theorem [8] we arrive at the inversion formulas 
Corollary 14. The solution to the inverse problem is given by: 



e^^k' -. 




k 

>>/ = 
j=0 


2 


C-l 
Cl 

Ck 


Co 

C2 
Cfe+1 


Cl 
C3 

Cfc+2 ■ 


• Cfe_i 

• Cfc+l 

C2fc 


, k-- 


= 0,1,.. 


. ,n- 


-1, 


A^i 



Proof. By Theorem [8] e^*' = -P2fc[0]. From Theorem [13] we obtain 

Q2k{x) 



P2k[0] 



-dfi{x), 



9 



which, in view of the determinantal formula presented in Theorem [T3l imphes the 
first formula. Likewise, to get the second formula we note 

and use Theorems 1131 and [8l to justify the claim. 

D 

Remark 15. Although we are not going to give a complete argument here, it is easy 
to show using for example the technique of Lemma 5.4 in [T] that all determinants 
appearing in the inversion are strictly positive. Also, it is worth noticing that all 
these determinants are in fact minors of the moment matrix /. 

Example 16. Here we solve the inverse problem for n = 1 and n = 2. If n = 1, by 
corollarv 1141 we get 

e^- - ^ 



1 -c_i. 

^0 



So 



and 



If n=2, we have 



xi = Jlogc_i = ^log(^) (2.29) 

2c_i 2 

mi = = — . (2.30) 

Co Ai 



e 



A} ci 



Therefore, 



g2x2 ^ ^ ^ 
^0 



1 c-ici-cl 1 aia2{Xi-X2)' 

xi = - log = - log ■ 



Also, we have 



2 ci 2 AiA2(aiAi + a2A2)' 

1 1 aiA2 + a2Xi 

"^ = 2^°S^"^ = 2^°S A1A2 • 



2 2(c_iC2 - CqCi, 

1711+1712 = — (t(c_iC2 - CoCl) = 5 

A^ C0C2 - Cf 

2c_i 2(aiA2 + a2Ai) 

1712 



(2.31) 



(2.32) 



Co AiA2(ai + 02) ' 



Therefore, 



_ 2ci(c_ici -eg) _ 2(aiAi +a2A2) , . 

co(coC2 - cf ) AiA2(Oi + 02) 

These formulas (after shifting Aj -^ —Xj, aj — > '^djXj ) are identical to what is given 
on page 246 in [l]. 
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